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^ ■ Abstract 
< 

Heavy quarks are copiously produced in ultra-peripheral heavy ion 

. collisions. In the strong electromagnetic fields, cc and bb are produced 

I by photonuclear and two-photon interactions. Hadroproduction can also 

I occur in grazing interactions. We calculate the total cross sections and the 

\^ ■ quark transverse momentum and rapidity distributions, as well as the QQ 

' invariant mass spectra from the three production channels. We consider 

. AA and pA collisions at the Relativistic Heavy Ion Collider and the Large 
Hadron Collider. We discuss techniques for separating the three processes 

Qh! and describe how the AA to pA production ratios might be measured 

p ' accurately enough to study nuclear shadowing. 
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I. INTRODUCTION 



In ultra-peripheral heavy ion collisions, heavy quarks can be produced in electro- 
magnetic or hadronic interactions. Electromagnetic production occurs through strong 
electromagnetic fields which can interact with a target nucleus in the opposing beam 
(photoproduction) or with the electromagnetic field of the opposing beam (two-photon 
reactions) and produce hadronic final states, including heavy quark pairs. Hadroproduc- 
tion of heavy quark pairs is also possible in grazing interactions. Since photon emission 
is coherent over the entire nucleus and because the photon is colorless, the three channels 
can be distinguished by the presence of zero, one or two rapidity gaps in the events and 
by whether or not the nuclei dissociate. 

Many types of ultra-peripheral collisions have been considered [0. Photonuclear 
interaction studies have included Coulomb dissociation and coherent vector meson 
production Final states studied in two-photon interactions have included produc- 
tion of lepton pairs [H , single mesons and meson pairs M as well as production of the 
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Higgs boson and other exotica 0]. Although the hst of experimentally observed channels 
is currently short, as RHIC gears up new results should come quickly and measurements 
of heavy quark production in ultra-peripheral collisions may not be too far off. 



We build on previous calculations of heavy quark photoproduction P|-pD| and two- 
photon production |Tl| in heavy ion collisions. For the first time, we consider resolved 
photon processes. We use modern parton distribution functions along with up-to-date 
accelerator species and luminosities. We also consider hadroproduction in grazing colli- 
sions and compare the three channels. Finally, since useful measurements of shadowing 
will require high accuracy, we consider the uncertainties inherent in these calculations 
and methods to control them. 



We will compare the total cross sections, heavy quark transverse momentum, px, and 
rapidity, y, distributions, and the QQ pair invariant mass, M, distributions in all three 
channels. We also discuss pA collisions to see how a comparison of photoproduction 
in AA (7^4) and pA (effectively 7p) can be used to study nuclear effects on the parton 
distribution functions (shadowing). 



For consistency, we will use parallel approaches to all three calculations with the 
same quark mass, parton distribution and QCD scale. For the charm quark mass, rric, 
this is somewhat problematic because studies of different production channels seem to 
prefer different values. Hadroproduction calculations have typically best fit the data 



with a relatively light rric, ~ 1.2 — 1.4 GeV |T^, while photoproduction studies have 



generally favored higher values, 1.5 — 1.8 GeV [13,14]. The limited two-photon data 
favors an intermediate value, around 1.6 GeV |TB|,[TB|. For the bottom quark, the typical 
mass range considered is 4.5 — 5.0 GeV. Hadroproduction calculations with mi, = 4.75 
GeV underpredict the bb cross section observed in 1.8 TeV pp collisions [|l3], suggesting 
that a smaller mass might be preferred [jl^ . Some exotic b production mechanisms have 
been suggested to explain the excess |jl8|,|l9[. Lattice studies suggest lower c and b quark 



masses, 1.1 — 1.4 GeV for charm and 4.1 — 4.4 GeV for bottom (the pole masses are 
somewhat higher) [^. We use nic = 1.2 GeV and rrif, = 4.75 GeV throughout this 
paper but we will also discuss the effect of varying the mass and scale. The QCD scale 
entering the running coupling constant, asiQ"^), and the parton distribution functions 
are proportional to the transverse masses, - 



4m^ for charm and m|. for bottom |12]. 



Table H gives recent estimates of nucleon-nucleon center of mass energies, y/S, and 
luminosities for AA and pA collisions at RHIC E^ and the LHC MM. The LHC 



luminosities assume that two experiments take data and that there is a 125 ns bunch 
spacing. At RHIC, the pA energies are the same as the AA energies while the luminosities 
are taken to be the geometric mean of the AA and pp luminosities {Cpp = 1.4 x 10^^ 
cm~^ s~^)- There is also the possibility of pA collisions at the LHC. However, at the 
LHC, the proton and the ion must have the same magnetic rigidity. Because protons 
and ions have different charge to mass ratios, they have different per nucleon energies 
and the center of mass is no longer at rest in the lab. These pA collisions are then 
at somewhat higher per nucleon energies than the corresponding AA collisions. At the 
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LHC, the pA collision rates can exceed 200 kHz although some experiments may need 
to run at a lower luminosity. Because RHIC is a dedicated heavy ion accelerator, it is 
likely to run a wider variety of beams than the LHC. At RHIC, dA collisions may be an 
alternative or supplement to pA. Except for the different initial isospin, most of the pA 
discussion should also hold for dA because shadowing is small in deuterium. 



Section |T| discusses photoproduction of heavy quarks in AA and pA collisions. Sec- 
tion m covers hadroproduction in peripheral AA and minimum bias pA collisions. Sec- 
tion |V] considers 77 QQ production. Section |V] compares our results for the three 
channels while section |V| is dedicated to a discussion of how to disentangle the produc- 
tion channels experimentally. In section |VII|, we draw our conclusions. 



II. PHOTOPRODUCTION 



Photoproduction of heavy quarks occurs when a photon emitted from one nucleus 
fuses with a gluon from the other nucleus, forming a QQ pair [p|,|9|,p^ ( "direct" produc- 
tion), as in Fig. |l|(a). The photon can also fluctuate into a state with multiple qq pairs 



and gluons, i.e. \n{qq)m{g)) . One of these photon components can interact with a quark 
or gluon from the target nucleus ("resolved" production), as in Figs. |I|(b)-(d) The 
photon components are described by parton densities similar to those used for protons 
except that no useful momentum sum rule applies to the photon 



At leading order (LO), the partonic cross section of the direct contribution is pro- 
portional to aas{Q'^)eQ, where as(Q^) is the strong coupling constant, a = e^/hc is the 
electromagnetic coupling constant, and eg is the quark charge, Cc = 2/3 and e?, = —1/3. 
The resolved partonic cross section is proportional to a'^{Q^). Even though the resolved 
partonic cross sections are larger than the direct partonic cross section, the smaller flux 
of quarks and gluons from the photon suggests that the resolved contribution should be 
smaller than the direct component. 



The cross sections are calculated using the Weizsacker- Williams virtual photon flux, 
modern parameterizations of the target gluon and quark distributions and the LO par- 
tonic cross sections. Newer parton distributions are considerably softer than the fiat, 
scaling parameterizations used earlier We require that the photoproduction not 

be accompanied by hadronic interactions 0]^. This could be done by restricting the 
impact parameter, b, to greater than twice the nuclear radius, Ra- Here we weight the 
6-dependent photoproduction probability by the 6-dependent hadronic non-interaction 
probability. 



Direct QQ pairs are produced in the reaction 'j{k) + N{P2) Qipi) + Q{P2) + ^ 
where k is the four momentum of the photon emitted from the virtual photon field of 
the projectile nucleus, P2 is the four momentum of the interacting nucleon in ion A, 
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and pi and p2 are the four momenta of the produced Q and Q. The photons are almost 
real. Their slight virtuality, Iq'^I < {Tic/RAf', is neglected. 

On the parton level, the photon-gluon fusion reaction is 7(fc) + g{x2P2) — > Q{pi) + 
Q{P2) where X2 is the fraction of the target momentum carried by the gluon. The LO 
QQ photopro duct ion cross section for quarks with mass mq is 



(Pa 

, ]^ = 7ias{Q'^)aeQBQED{s, h, ui)5{s + h + ui) (1) 



where 
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Here as{Q^) is evaluated to one loop at scale Q^. The partonic invariants, s, ti, and 
ui, are defined as s = {k + X2P2Y, ti = {k — piY — mq = {X2P2 — P2Y — ^q, and 
ui = {X2P2 — PiY — = {k — P2Y — niQ. In this case, s = 4:kjLX2mp where •jl is 
the Lorentz boost of a single beam and rrip is the proton mass. Since k ranges over 
a continuum of energies up to -Ebcam = iL^p, we define Xi = k/Pi analogous to the 
parton momentum fraction where Pi is the nucleon four momentum. For a detected 
quark in a nucleon-nucleon collision, the hadronic invariants are then S = (Pi + P2Y, 
Ti = (P2 - pif - m^, and U, = (Pi - p,f - m^. 

We label the quark rapidity as yi and the antiquark rapidity as y2- The quark 
rapidity is related to the invariant Ti by Ti = —\/SmTe~'^^ where iriT = \Jpt + rnq. 

The invariant mass of the pair can be determined if both the Q and Q are detected. The 
square of the invariant mass, = s = 2m^(l + cosh(?/i — j/2)), is the partonic center of 
mass energy squared. The minimum photon momentum necessary to produce a QQ pair 
is fcmin = M^/A-fLmp. At LO, X2 = {mT/VS){ey^+ey^) and xi = {mT/VS){e-y' +6-^^). 
We calculate xi and X2 as in an A^A^ collision and then determine the flux in the lab frame 
for k = Xi7imp, equivalent to the center of mass frame in a collider. The photon flux 
is exponentially suppressed for k > '^iJic/Ra-, corresponding to a momentum fraction 
Xi > hc/nipRA- The maximum 7A^ center of mass energy, JSjn, is much lower than the 



hadronic v^- Note that S^n = W-y^, the typical notation for HERA. For consistency, 
we use "s/S* notation for all three processes. 

The cross section for direct photon-nucleon heavy quark photoproduction is obtained 
by convoluting Eq. (|I]) with the photon flux and the gluon distribution in the nucleus 
and integrating over k and X2, 

,2 ^XW _2 r dk^ t '^Ft{x2.Q\h^z)s^^ , (3) 



dTidUi(Ph J Jk^in dkd?h Jx2^^^ X2 ^ 111 dtidui 
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where d^N^/dkd% is the differential photon flux from one nucleus (our final results will be 
integrated over h > 2Ra) and z is the longitudinal distance. The factor of two in Eq. (||) 
arises because both nuclei emit photons and thus serve as targets. The incoherence 
of heavy quark production eliminates interference between the two production sources 
|28[] . Four- momentum conservation gives a;2„ji„ = —Ui/{S + Ti) in terms of the nucleon- 
nucleon invariants. The equivalent hadronic invariants can be defined for photon four 
momentum k as S^n = {k+P2y, Ti^-yN = (-P2— Pi)^— '^g, and ?7i,^7v = {k—pi)'^—m^ |2^. 
The partonic and equivalent hadronic invariants for fixed k are related by s = X2S^n, 
ti = Ui^^N, and Ui = X2Ti^'yN- 



We now turn to the resolved (hadronic) contribution to the photoproduction cross 
section. The hadronic reaction, QQX, is unchanged, but now, prior to the inter- 

action with the nucleon, the photon splits into a color singlet state with some number 
of qq pairs and gluons. There are a few photon parton distributions available [30-3|]. 
None of them can be definitively ruled out by the existing data on the photon structure 
function p5|j36[] . As expected, F^{x,Q'^) = F^{x,Q'^) flavor by flavor because there are 
no "valence" quarks in the photon. The gluon distribution in the photon is less well 
known. We use the GRV LO set [^]. Its gluon distribution is similar to most of the 
other available sets |13T1 , |33| , |3^ . Only the LACl set has a higher low- a; gluon density. 



up to an order of magnitude larger than the others. 



On the parton level, the resolved LO reactions are g{xk) + g{x2P2) Q{Pi) + 
Q{p2) and q{xk) + g(x2-P2) Q{Pi) + Q{P2) where x is the fraction of the photon 
momentum carried by the parton. The LO diagrams for resolved photoproduction, 
shown in Fig. [I|(b)-(d), are the same as for hadroproduction except that one parton 



source is a photon rather than a nucleon. The LO partonic cross sections are |37 
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dtidui 
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-Bqed{s, ii, ui 




2tiUi 



S{s + ti + ill) , 



(4) 
(5) 



where s = {xk+X2P2y,ii = (x/c— pi)^— mg, andwi = {x2P2—pi)'^—mQ. The (7(7 partonic 
cross section, Eq. (^, is proportional to the photon-gluon fusion cross section, Eq. (|T]), 
with an additional factor for the non-Abelian three-gluon vertex. The qq annihilation 
cross section has a different structure because it is an s-channel process with gluon 
exchange between the qq and QQ vertices. The gg reactions are shown in Figs. |l](b)-(c); 
Fig. Kc) is the non-Abelian contribution. The diagram is shown in Fig. |l](d). Modulo 
the additional factor in the gg cross section, the resolved partonic photoproduction cross 
sections are a factor as{Q^)/aeQ larger than the direct, •yg, partonic photoproduction 
cross sections. Despite this, the resolved component is still smaller than the direct 
component. 

The cross section for resolved photoproduction is 
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where fcmin is defined as before. Since k is typically larger in resolved than direct photo- 
production, the average photon flux is lower in the resolved contribution. 

The nuclear parton densities F/^{x, Q^, b, z) in Eqs. (||) and ® can be factorized 
into X and independent nuclear density distributions, position and nuclear-number 
independent nucleon parton densities, and a shadowing function x, Q^, -z) that 

describes the modification of the nuclear parton distributions in position and momentum 



space. Then |8[f| 



(7) 



where fl^ {x,Q ) is the parton density in the nucleon. We evaluate the MRST LO 
parton distributions [43| at = a^m^ where a = 2 for charm and 1 for bottom. 



In the absence of nuclear modifications, S^{A,x,Q'^ ,b, z) = 1. The nuclear density 
distribution, pA^b, 2;), is a Woods-Saxon shape with parameters determined from electron 
scattering data [Q]. Although most models of shadowing predict a dependence on the 
parton position in the nucleus, in this photoproduction calculation we neglect any impact 
parameter dependence. Then the position dependence drops out of 5*. We employ the 
EKS98 shadowing parameterization available in PDFLIB PH| for S** 7^ 1. 



The full photoproduction cross section is the sum of the direct and resolved contri- 
butions |p!3| . 



d'^n - 
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dTidUid^b 
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The total cross section is the integral over impact parameter and the hadronic invariants 
Ti and Ui, 



■^A^QQX 



dTi dUi d% 



dV 



jA^QQX 

dTidUid^b 



(9) 



When 5* = 1, the impact-parameter integrated cross section in Eq. 
Including shadowing makes the dependence on A nonlinear. 



scales with A. 
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The photon flux is given by the Weizsacker-Wilhams method. The flux from a charge 
Z nucleus a distance r away is 



dkd^r n'^kr'^ 



k!{w) + - 

1l 



(10) 



where w = kr/'ji and Ko{w) and Ki{w) are modified Bessel functions. The photon flux 
decreases exponentially above a cutoff energy determined by the size of the nucleus. In 
the lab frame, the cutoff is fcmax ~ 1l^c/Ra- In the rest frame of the target nucleus, the 
cutoff is boosted to -Emax = (27I, — l)hc/RA- Table y shows the beam energies, -Ebeam, 
Lorentz factors, 7^, fcmax, and -Emax, as well as the corresponding maximum center of 
mass energy, ^/Smal = \/2-E'max^p, for single photon interactions with protons, 7p — ^ QQ 



46[| . At the LHC, the energies are high enough for tt photoproduction |jTO| 



The total photon flux striking the target nucleus is the integral of Eq. (|T0|) over the 
transverse area of the target for all impact parameters subject to the constraint that 
the two nuclei do not interact hadronically 0]. This must be calculated numerically. 
However, a reasonable analytic approximation for AA collisions is given by the photon 
flux integrated over radii r > 2i?^. The analytic photon flux is 



dN^ 
dk 



2Z'a 
nk 



'111 



where w^"^ = 2kRAllL- We use the more accurate numerical calculations here. The 
difference between the numerical and analytic expressions is typically less than 15%, 
except for photon energies near the cutoff. The analytical and numerical photon fluxes 
differ most for hb production at RHIC. 

The photoproduction distributions are shown in Figs. |]-| for the largest nuclei at 
each energy, gold for RHIC and lead for LHC. Since y/S^^ is close to the bb production 
threshold for iodine and gold beams at RHIC, the pt and mass distributions for these 
nuclei are narrower than those with oxygen and silicon beams. The direct photopro- 
duction results are reduced by a factor of two on the figures to separate them from the 
total. There are two curves for each contribution, one without shadowing, S** = 1, and 
one with homogeneous nuclear shadowing, 5* = EKS98. When the effects of shadowing 
are small, the curves are indistinguishable. 

Shadowing has the largest effect on the rapidity distributions in Figs. i|(b)-|^(b). In 
these calculations, the photon is emitted from the nucleus coming from positive rapidity. 
Then j/i < corresponds to < 7LX2Wip in the center of mass (lab) frame. If the 
photon emitter and target nucleus are interchanged, the resulting unshadowed rapidity 
distribution, S"* = 1, is the mirror image of these distributions around yi = 0. The Q 
and Q distributions are asymmetric around yi = 0. The resolved contribution is largest 
at rapidities where the photon momentum is small. The resolved rapidity distributions 
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are narrower and shifted to larger negative rapidity than the direct contribution. The 
average mass and transverse momentum for the resolved component are smaller than 
for the direct ones. The total heavy quark rapidity distributions are the sum of the 
displayed results with their mirror images when both nuclei emit photons. This factor 
of two is included in the transverse momentum and invariant mass distributions. 



The impact-parameter integrated total direct and resolved photoproduction cross 
sections are given in Table |T| for charm and Table |^ for bottom. The difference 
between the 5* = 1 and = EKS98 calculations is due to shadowing. The change 
is about 10% for cc with gold at RHIC, rising to 20% for cc with lead at the LHC 
For bb, the shadowing effect is smaller, about 5%. By selecting charm production in a 
limited rapidity and px range, it would be possible to enhance the effect of shadowing 
slightly but the 10% and 20% effects at RHIC and the LHC are useful benchmarks of 
the accuracy needed for a meaningful measurement. 



The cross sections vary by orders of magnitude between the lightest and heaviest 
targets, primarily due to the changing Z^. At RHIC, (Zau/^o)^ is almost exactly the 
ratio of the bb production cross sections. For bb production at the LHC, the lead to 
oxygen cross section ratio is ~ 60, less than {Zpy,/ZoY- 

After adjustment for different parton distributions, quark masses, beam energies and 
species as well as the resolved contributions, the photoproduction results are comparable 
with previous studies. The direct cc cross sections in Table fT^ are almost identical to 
those found by Baron and Baur despite significant differences in gluon structure 
function and quark masses {e.g. rric = 1.74 GeV in Ref. rather than rric = 1.2 GeV). 
The quark mass difference can be compensated by the newer parton distributions since 
the low-x MRST gluon density is much higher than that of the older Duke-Owens parton 
distributions used by Baron and Baur. Our bb cross section is about three times higher 
than Greiner et al. at the same invariant mass although they use a smaller b quark 

mass, 4.5 GeV, due to the larger low-x MRST gluon distribution. 

On the other hand, nonperturbative QCD calculations of direct photoproduction can 
yield very different rates. A colored glass model predicts QQ production cross sections 
in Pb-|-Pb collisions at the LHC of 800 mb for cc and 100 mb for bb about half 
the perturbative cc cross section but 140 times the hb cross section. The ratio of charm 
to bottom cross sections in this formulation is extraordinarily small compared to the 
perturbative results in Tables ^ and 0. 

We now discuss the resolved rates. At RHIC the cc and hb resolved contributions are 
~ 2 and 6% of the total. The contribution is small because the available 7p energy is 
relatively close to threshold: the resolved center of mass energy is smaller than i/Smax- 
At the LHC, the resolved contributions are ~ 15 and 20% of the total charm and 
bottom photoproduction cross sections respectively, comparable to the shadowing effect. 
Interestingly, the resolved qq component is considerably greater than qq annihilation in 
hadroproduction because, like the valence quark distributions of the proton, the photon 
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quark and antiquark distributions peak at large x. However, the peak of the photon 
quark distribution is at higher x than in the proton. These large photon momentum 
fractions increase the qq contribution near threshold. The total cc resolved cross section 
at RHIC is (35 — 50)% qq while the hb resolved contribution is (80 — 90)% qq. The 
relative resolved contribution remains larger than qq annihilation in hadroproduction at 
the LHC where the qq contribution to the hadroproduction cross section is 4% for charm 
and 10% for bottom. The LACl gluon distribution predicts a considerably larger 
resolved contribution with a smaller relative qq contribution. 



At fixed target energies, near threshold, the resolved component of charm photo- 
production is relatively unimportant ||5l|. However, at higher 7p energies such as those 
at HERA, the resolved component becomes important. RHIC is an intermediate case, 
with an average 7p energy less than 1/5 of that at HERA. We find that, at RHIC, the 
resolved component is a (2 — 6)% effect. At the LHC, where the energies are gener- 
ally higher than at HERA, it is (15 — 20)% of the total cross section, concentrated at 
large rapidities. As was previously mentioned, much early photoproduction data favors 



heavier quark masses, 1.5 — 1.8 GeV, although newer HERA results are ambiguous ||52 |. 
Because we use a lower m^, our cross sections are higher than a direct comparison to 
HERA would indicate. Changing iric affects the overall cross section and the quark pt 



spectrum for < 5 GeV [13 



As is typically done, we include all QQ pairs in the total cross sections and rates even 
though some of these pairs have masses below the HH threshold where HH = DD, BB. 
Photoproduction is an inclusive process; accompanying particles can combine with the 
Q and Q, allowing the pairs with M < 2m h to hadronize. We assume the hadroniza- 
tion process does not affect the rate. Including all pairs in the total cross section is 
presumably an even safer assumption for hadroproduction because there are more ac- 
companying particles. On the other hand, 77 interactions should have no additional 
particles since the interaction is purely electromagnetic and occurs away from the nuclei 
in free space. Section ^ will discuss this in more detail. 



One way to avoid some uncertainties due to higher order corrections is to measure 
shadowing by comparing the pA and AA photoproduction cross sections at equal photon 
energies since the parameter dependence cancels in the ratio a{AA)/a{pA). In the equal 
speed system, equal photon energies correspond to the same final-state rapidities. In pA 
collisions, the photon almost always comes from the nucleus due to its stronger field. 
Thus the pA rates depend on the free proton gluon distribution. The photon fluxes are 
different for pA and A A because the minimum radii used to determine ujji are different: 
2Ra in AA compared to Ra + Tp in pA where Vp is the proton radius. There are a 
number of ways to define the proton radius. We use the hadronic radius, Vp ^ 0.6 fm, 
determined from photoproduction data As we will show, our results are not very 
sensitive to Vp. 

In a detailed calculation, the hadronic interaction probability near the minimum 
radius depends on the matter distribution in the nucleus. Our calculations use Woods- 
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Saxon distributions with parameters fit to electron scattering data. This data is quite 
accurate. However, electron scattering is only sensitive to the charge distribution in the 
nucleus. Recent measurements indicate that the neutron and proton distributions differ 
in nuclei ||5^. This uncertainty in the matter distribution is likely to limit the photon 



flux determination. 



The uncertainty in the photon flux can be reduced by calibrating it with other 
measurements. Further studies of the matter distributions in nuclei and the total 
ion-ion cross sections, as well as comparisons with known photoproduction and two- 
photon processes, can help pin down the photon flux. For example, the cross section 
for 7p — > pp is known to 10% from measurements at HERA [Q. The 7^4 — > VA 
data are limited to energies lower than but still comparable to those available at RHIC. 
A combination of lower energy RHIC data and a judicious extrapolation could fix the 
calibration. Studies of well known two-photon processes, like lepton production, can 
also help refine the determination of the photon flux. With these checks, it should be 
possible to understand the photon flux in pA relative to AA to better than 10%, good 
enough for a useful shadowing measurement. 

Our pA results are calculated in the equal-speed frame. This is appropriate for RHIC 
where \fS is the same in pA and AA interactions. At the LHC, the proton and nuclear 
beams must have the same magnetic rigidity and, hence, different velocities and per- 
nucleon energies. At the LHC, the equal-speed frame is moving in the laboratory frame 
so that the maximum pA energy is larger than the A A energy. The 71, and given 
for pA at the LHC in Tables H and |I| are those of the equal-speed system. The pA total 
cross sections for QQ production are given in Table [V|. 

It is easy to compare pA and AA results at RHIC. For = 1, the only difference is 
in the impact parameter range: b > 2Ra in AA compared to 6 > Ra + Tp in pA. We 
compare the numerical AA results presented in Tables |T| and |V| with those using the 
analytical photon flux for pA in Table |V|. To normalize the photon cross sections, we 
divide the AA cross section by 2A because there is only a single photon source in pA and 
the proton target is a single nucleon. Due to the reduced minimum impact parameters, 
a{AA)/2A is (22 - 37)% lower for cc and (37 - 65)% lower for bb with rp = 0.6 fm. 
The results differ least for small A where Ra — Vp is reduced. Increasing Vp by 50% 
would decrease the cc cross sections (6 — 3.5)% and the bb cross sections by (11 — 7)%. 
Changing Vp has the largest effect for small A since Ra/^p is smaller. The hb differences 
are largest because the energy is close to production threshold and high k corresponds 
to small impact parameter. 

At the LHC, the maximum pA and AA energies differ, adding another variable to 
the comparison. Here, we compare a{AA)/2A directly with (y{pA), using the analytic 
flux for both to maximize the parallelism. For equal pA and AA energies, the different b 
range is less significant at the LHC. The decrease in AA relative to pA is (18 — 23)% for 
cc and (24 — 30)% for bb. As at RHIC, the effects are larger for heavier nuclei. Increasing 
Tp by 50% is only a (2 — 5)% effect with the bigger change for bb. 
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The effect of the higher pA energy is significant. At the lower A A beam energy, 
the pA cross section is reduced (23 — 32)% for cc and (30 — 40)% for bb relative to the 
maximum pA energy. The major difference between the cc and bb results is likely due to 
the smaller x values probed in cc production. The increased energy has a larger effect 
on the cross section than the change in the minimum impact parameter. The energy 
dependence illustrates the desirability of obtaining energy excitation functions in both 
pA and AA interactions. 

Tables |V^ and |V11| give the total monthly AA and pA QQ production rates assuming 
a 10^ s/month run. The cc rates are quite high. Of course, hadronization, branching 
ratios and detector acceptances will all reduce the observed rates. The bb rates are only 
significant at LHC However, RHIC is a dedicated heavy ion facility, originally planned 
for more than 30 weeks of operation per year compared to 4 weeks of heavy ion running 
at the LHC and should thus accumulate more data than the tables indicate. 



Hadroproduction of heavy quarks in heavy ion collisions has been considered by 
many authors, see e.g. P^ , pD|J^ , |5D[ ] and references therein. Here, we consider the special 
characteristics of heavy quark production in peripheral heavy ion collisions. Two aspects 
meriting our attention are the form of the parton densities at the nuclear surface and the 
overall appearance of the collision — could hadroproduction mimic a photon-nucleus or 
two-photon interaction? This section addresses the first issue while section ^ considers 
the second. 

At leading order, heavy quarks are produced via the reaction N{Pi) + N{P2) 
Q{pi) + Q{p2) + X. The LO parton reactions are g{xiPi) + g{x2P2) — ^ Q{Pi) + QiP2) 
and g(xiPi) + g(x2-P2) Qipi) + Q{P2)- The LO partonic cross sections are those 
of Eqs. (ID and (|) with s, ti, and Ui replaced by s, ti and Ui. Thus the partonic 
couplings are a factor as{Q^)/(yeQ larger than direct photoproduction but have the same 
coupling strength as resolved photoproduction, as discussed in the previous section. The 
partonic invariants, s, ti, and Ui, are now s = {xiPi +X2P2)^, ti = {xiPi —piY — f^Q = 
{X2P2 - P2y - "^Q, and Ml = {X2P2 - Pi)^ - = {xiPi - ^2)^ - rrig. 

The hadronic heavy quark production cross section is the convolution of two nuclear 
parton distributions with the partonic cross sections. 



III. HADROPRODUCTION 



s- 



,2 "AA-^QQX 



I dzdz'^^ F/(xl,Q^r,^)F/(x2,QM6-r1,^')^^ 
J Xi X2 \ " ^ 

+ E {F^{x^,Q\T,z)Ff{x2,Q\\b-r\z') 




Q ^gg 
dtidu] 



dTidUidVr 
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+ F^^{x„Q',f,z)Ff{x„QMb-flz')}s'-^) . (12) 



If the quark is detected, the hadronic invariants are again S = (Pi + ^2)^, ^1 = (-P2 — 
Pi)"^ — itlq, and Ui = (Pi — pi)"^ — ttlq [0. The partonic and hadronic invariants are 



now related by s = X1X2S, ti = xiUi, and ui = X2T1. Four-momentum conservation 
at leading order gives X2^i„ = —xiUi/ {xiS + Ti). We again perform a fully leading 
order calculation using the MRST LO parton distributions and a one-loop evaluation of 

Since both the projectile and target partons come from nuclei, we consider the spatial 
dependence of the nuclear parton distributions, as suggested by shadowing measurements 
|55[] , in more detail. Shadowing in peripheral collisions could then be significantly dif- 



ferent than in central collisions. We consider three scenarios: no shadowing, = 1; 
homogeneous shadowing with the EKS98 parameterization, S** = EKS98; and inhomo- 
geneous shadowing, 5"* = EKS986, with spatial dependence included, as described in 
Refs. 



To obtain the AA cross sections, we integrate Eq. (|T2D over d'^r and d'^b. When 
5** = 1 or EKS98, the only b dependence is in the nuclear density distributions and 

aAAOCdNN J d%d'^rTA{r)TA{\b-f\) (13) 



where T^(r) = / dzpA^r, z) is the nuclear profile function and aj^^ is the QQ production 
cross section in an A^A^ collision. When Eq. ([l^) is integrated over all b, (Jaa oc ctatatA^ 
for 5** = 1. Homogeneous shadowing makes the A dependence nonlinear so that the 
integrated cross section is effectively aAA oc a^NA'^"' where a can be determined from 
the A dependence of the EKS98 parameterization. Of course when S** = EKS986, the 
full integral over b and r is needed in Eq. (|1^). In this section the impact parameter 
integral is restricted to 6 > 2P^ to be consistent with the photoproduction results. 

This calculation treats the nuclei as a continuous fluid, neglecting the lumpiness 
due to the individual nucleons. When b > 2Ra, only a handful of nucleon-nucleon 
collisions can occur, resulting in significant statistical fluctuations. These fluctuations, 
although unimportant for heavy quark production, may affect some of the observables 
used for impact parameter determination such as transverse energy or charged particle 
multiplicity. 

The hadroproduction distributions for the heaviest ions are shown in Figs. P-^. There 
are three curves for each colliding system. When the small x region is probed, the S"* = 1 
result is the highest, EKS98 is the lowest, and EKS986, with b > 2Ra, is between the 
other two. If we consider minimum impact parameters much larger than 2P^, the 
EKS986 curves would move even closer to the = 1 results. The order is reversed for 
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bb production at RHIC because x ~ M/y/S ~ 0.05 - 0.1 for 2mb < M < 20 GeV, in the 
gluon antishadowing region of the EKS98 parameterization. Thus the shadowed results 
he above those with 5* = 1. The quark pt and QQ pair mass distributions are harder 
than the photoproduction results in Figs. since now a/S* is the full nucleon-nucleon 
center of mass energy. In hadroproduction, the quark rapidity distribution is symmetric 
around yi = 0. 



The cc AA cross sections and rates are given in Tables [VII1| and The bb cross 
sections and rates are in Tables ^ and ^ respectively. Shadowing has less effect on the 
total cross sections than in our previous calculations [p8|,p9| which used earher shadowing 



parameterizations [l56| , |57| with stronger gluon shadowing at low x and weaker gluon 
antishadowing. 



Shadowing effects depend on the final state rapidity and pair mass. The region away 



from yi = tends to be most sensitive to the shadowing parameterization At 
RHIC, the effect grows with rapidity because at ?/i ~ a; is not small and shadowing is 
weak (charm) or x is in the antishadowing region (bottom). Higher positive rapidities 
correspond to lower X2 for the target (stronger low x shadowing) and larger xi for the 
projectile (the EMC region), increasing the effect. 

For bb production, the x region moves from antishadowing at ?/i = to shadowing as 
yi increases. In Au+Au colhsions at RHIC, a{S' = EKS98) /a{S' = 1) = 1.22 at ?/i = 
and 0.917 at yi = 2.5. At the LHC, the cross section varies less with rapidity because 
both the target and projectile momentum fractions are in the shadowing region. In fact 
shadowing tends to decrease slightly with rapidity because 5"' increases more with xi 
than it decreases with X2. 



Nuclear shadowing is more important at the LHC than at RHIC because the higher 
energy collisions probe smaller x values. Shadowing is also larger for cc than bb because 
the lighter charm quark is produced by lower x partons than b quarks. Homogeneous 
shadowing reduces the cc total cross section by (3 — 5)% at RHIC while at the LHC the 
reduction is (18 — 34)%. As expected, the reduction is largest for the heaviest nuclei. 
With the inhomogeneous shadowing parameterization, EKS986, the cc cross section is 
reduced by (1 — 2)% at RHIC and ~ 10% at LHC. The bb cross section is antishadowed 
at RHIC by (6-20)% using EKS98, reduced to (3-8)% with EKS986. At the LHC, the 
bb cross section is reduced by (10 — 19)% with EKS98 and ~ 6% with EKS986. At large 
b, any inhomogeneity reduces the effect of shadowing. The larger the impact parameter 
cut, the closer the peripheral AA events mimic hadroproduction in free proton collisions. 

A direct comparison with the photoproduction cross sections in Tables |T| and fV| 
shows that the only rate comparable to hadroproduction at RHIC is cc photoproduction 
in Au+Au collisions when the same impact parameter range is considered. The photo- 
production cross section is ~ 25 mb, comparable to the 33 mb hadroproduction cross 
section. With lead beams at the LHC, the cc photoproduction cross section is actually 
larger than the hadroproduction cross section. The bb photoproduction cross section is 
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always much lower than that of hadropro duct ion. 



The hadropro duct ion cross sections tend to be larger than photoproduction cross 
sections for several reasons. In hadroproduction, the full is available for heavy 
quark production while coherent photon emission requires ^S^n ^ V^- The lower 
energy reduces the cross section and restricts the x range of the gluons taking part in 
the interactions. Thus hadroproduction probes smaller x values than photoproduction. 
At low X the gluon densities are larger than the photon flux. 



The minimum bias pA results with = 1 and = EKS98 are presented in Ta- 
bles pCIH and pCIIl] for charm and bottom respectively. We report the minimum bias pA 
cross sections only since it is difficult to select peripheral pA events. Shadowing is less 
important than in AA collisions since the pA cross section is linear in S"* while the AA 
cross section is quadratic in S^. The minimum bias cross section is proportional to A 
for = 1. A comparison of the RHIC minimum bias pA and the peripheral AA cross 
sections shows that the pA cross section without shadowing is equal to the AA cross 
section divided by (l/A) J d?hTAA(p > '^Ra), as expected. There is no corresponding 
factor of A for photoproduction so the hadroproduction pA cross sections are always 
bigger than the photoproduction cross sections in pA. Recent studies have shown that a 
comparison of hadroproduction in pA and pp collisions at the same energies can provide 
detailed information on nuclear shadowing effects |^|. In this case, there is no difference 



in flux between pA and AA collisions as there is in photoproduction. Such studies could 
provide important input to the AA collisions discussed here. 



The variations in the cross section due to quark mass and QCD scale are similar in 
both hadroproduction and photoproduction. However, the additional NLO corrections 
are larger in hadroproduction. Even at NLO the calculations do not always agree with 
data. The measured production rate in pp collisions at ^ = 1800 GeV is 2.9 ± 
0.2 ± 0.4 times the NLO calculation The reason for this discrepancy is unknown 
but some non-standard suggestions have been made p!8|JT9[1 . The discrepancy may also 
be due to an incomplete understanding of the hadronization process 



The major uncertainty for hadroproduction in peripheral collisions is the minimum 
impact parameter. There is no known method for effectively selecting very large-impact 
parameter hadronic events or, alternatively, collisions with a small but well defined 
number of participants. Zero degree calorimeters (ZDCs) can be used to select events 
with a small number of spectator neutrons but these come from a poorly defined range 
of impact parameters. For this reason, peripheral AA collisions may best be studied 
by comparing different processes in a ZDC-selected data set. These different processes, 
such as heavy quark and gauge boson production , could be used to compare nuclear 
parton densities for several species at a variety x and values. 
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IV. TWO-PHOTON PRODUCTION 



Heavy quark pairs can also be produced in purely electromagnetic photon-photon 
collisions. This process has been studied at e"'"e~ colliders. However, in ion colliders the 
cross sections are enhanced since the 77 luminosity increases as Z^. 



The 77 luminosity has been calculated by many authors |^ . Naively, it is given by 



the convolution of the photon fluxes from two ultrarelativistic nuclei. The photon flux 
from one nucleus is the integral of d?N^/ dkd^r in Eq. (0) over r excluding the nuclear 
interior. Not all the flux is usable because, when the nuclei actually collide, the two- 
photon interaction products will be lost amongst the much denser hadronic debris. The 
usable flux is limited by the requirement that the nuclei do not interact hadronically. 
This is typically done by requiring that h > 2Ra- However, when the photon energy is 
close to the kinematic limit, k^s.x ~ 1l^c/Ra, the flux is sensitive to the exact choice of 
Ra- To reduce the sensitivity to Ra, we calculate the probability, P{b), of a hadronic 
interaction as a function of impact parameter. 



P{b) = 1 — exp 



- J d'rTA{r)TB{\b - f\) 



(14) 



and use it to calculate the effective photon flux. Woods-Saxon density distributions 



3^ are used to calculate the nuclear overlap integral. The nucleon-nucleon total cross 
section, a]^^, is 52 mb at 200 GeV and 93 mb at 5.5 TeV We use the total 



cross section to exclude any interaction which could cause the nuclei to break up. This 
soft cutoff on the flux reduces the effective two-photon luminosity by a few percent for 
k -C ■jLfic/RA, rising to about 15% at the kinematic limit compared to a hard cutoff 
with Ra = 1.16(1 — 1.16A~'^^^)A^^^ fm. We also exclude the flux when the heavy quarks 
are produced inside one of the nuclei. Although these heavy quarks would probably 
survive, the resulting interactions are likely to disrupt the nucleus, making the collision 
appear hadronic. With these exclusions, the differential 77 luminosity is 

'"^'^ -^aa[_ / „ d%d^r4^^^Pi\b-f\). (15) 



dkidk2 Jb>RA Jr>RA dkid% dk2d'^r 



The two-photon center of mass energy, JS^^, is given by the photon energies, 5"^^ 



4:kik2- This 5*^^ is equivalent to W'^, commonly used in other two-photon calculations. 



The maximum yS^^ is 2/cinax = '^Il^/Rai given in Table This limit is a factor of 

(hc/nipRA)'^ smaller than v^, a factor of 10"'^ for gold or lead. Indeed, 2fcmax < 2mb 
for I+I and Au+Au collisions at RHIC. Thus heavy quark production in this channel 
should be much smaller than for photo- and hadropro duct ion. 

As in photoproduction, there are also direct and resolved contributions. Either one 



or both of the photons |15] can resolve itself into partons. At the parton level, the single- 
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resolved photon case is equivalent to photoproduction while the double-resolved photon 
situation is equivalent to hadropro duct ion. Both of these contributions are included 
here. The diagrams for all of these processes are shown in Fig. [10. 



The LO cross section for heavy quark production in direct two-photon interactions 
is also proportional to -Bqed, as in Eq. (|l]) for direct photoproduction, but with different 
couplings, 

(Pa 

s'^^—P- = 67ra^eQBQED{s,ti,Ui)S{s + ti + ui) , (16) 



dtidu 



where s = (fci + ^2)^ = S^^, ti = {ki — pi)^ — mg, and Ui = {k-i —pi)'^ — mQ. Here ti = Ui 
and Ml = Ti for a detected quark. The 77 QQ cross section is a factor of Qe^a/asiQ"^) 
smaller than the partonic 7(7 QQ cross section, Eq. ([1|). A factor of 3 comes from the 
three quark colors while another factor of 2 is from the spin averages. The ratio a/asiQ"^) 
reduces the cross section for two-photon production relative to photoproduction. The 
same two Compton diagrams apply to both two-photon production and photoproduction 
except that a second photon replaces the gluon in 77 interactions. 

The direct cross section is the convolution of the partonic two-photon cross section 
with the two-photon luminosity for photon energies ki and ^2, 

" dlidUi J dkidk2 dtidui 



The photon fluxes are exponentially suppressed for k > ■yih/RA- 

The resolved processes follow from the discussion in Section |T[ The cross section for 
singly resolved production of heavy quarks is 



^2^1-res _ rl rir rl"^ 

= 2 / dk.dk', I —j^f]i^2,Q')s'^ 

dlidUl J J X2 dkidk2 dtidui 



where the partonic invariants, s = {ki + X2k2Y = X2S^^, ii = (ki — piY — rriQ, and 
III = {x2k — pif — rriQ, are related to the total invariants by ii = Ui and ui = X2T1 
and k2 = X2k2. The partonic cross section, a^g, is the photoproduction cross section in 
Eq. (|I]). The cross section for double-resolved heavy quark production is 



_ rl rl rl r f rl"^ 

..-.QQ f fdX2 fax, dL f.^^^^Q2^f.^^^^Q2f±_ 

'^'^2 L dti 

;2 d'^a, 



q2 H-*QQ _ / AV.I AU 



dTidUi J J X2 J xi dk'idk'2 



+ E if^i^uQ')f^i^2,Q') + f^{xuQ')f^{x2,Q'))s 



QQ 



dtidui 



(19) 
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where ^2 = X2k2, k[ = xiki and the partonic invariants, s = {xiki + X2k2Y = xiX2S^j, 

ii = {xiki — pi)^ ""^Q) ui = {x2k2 —piY ~"^Q) related to the total invariants by 

ii = xiUi and ui = X2T1. The partonic cross sections, Oqq and agg, are given in Eqs. (^) 
and d). 

The full two-photon heavy quark cross section is the sum of all three contributions, 

^2 ^2 dir ciV^-'"''^^ f^'f^'""^7T 

o2 " ^ r'2 77->QQ I 02 77^QQ , 02 77^QQ /of)^ 



In Eqs. ( P^ ) and (|T^ the scale entering into the photon parton densities and Q!s(Q^) has 
been set equal to 4mQ due to the structure of the 77 luminosity. Changing the scale 
from 4mfn to 4mQ increases the single-resolved cross section by about 10%, while the 
double-resolved cross section changes by at most 2%. 



Figures pTUT^ show the corresponding cc and bb production distributions. The RHIC 
results are shown for Si+Si collisions since that is the largest A for which 2/cmax > 
2mb- The Pb+Pb results are shown for LHC The integrated cross sections for all the 
other nuclei are given in Tables pCIV| and pCV| . The quark pt and rapidity and the QQ 
pair invariant mass distributions are narrower for the heavier nuclei due to the lower 
associated ^Js^■ The production is mostly direct. The single and double resolved px 
and mass distributions are narrower than the direct results at all energies. 



The rapidity distributions are symmetric around yi = except for the singly resolved 
processes. Since either photon may be resolved, we give the single-resolved rapidity 
distributions in both cases. The total single-resolved rapidity distribution is the sum. 
We present both to be consistent with the direct and double-resolved calculations. The 
factor of two in the cross section, Eq. (pISD, is included in the single-resolved transverse 
momentum and invariant mass distributions. 



Direct 77 production dominates two-photon production of heavy quarks. For cc 
pairs, at RHIC the single-resolved cross section is (0.6 — 1.6)% of the direct production. 
Double-resolved production is (2.5 — 3.3)% of the single resolved. The single-resolved to 
direct ratio increases with A or, equivalently, decreasing \J S^^. Interestingly, the single- 
to double-resolved cross section ratios decrease with increasing A, presumably due to the 
larger qq component at lower energies, closer to production threshold where the quark 
distribution in the photon is dominant. At LHC, the single-resolved cross section is 4.5% 
of the direct production. The single-resolved component is higher at LHC because of the 
higher energy. However, double- resolved production is only 2.1% of the single resolved 
because the higher beam energy reduces x, increasing the photon gluon distribution but 
the gg contribution is too small to make up the difference. The situation changes for hb 
production. At RHIC, single-resolved production is (0.3 — 0.7)% of the direct component 
but double-resolved production is 20% of the single resolved. The lower energy strongly 
reduces the single-resolved cross section relative to the direct 77 but has less effect 
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on double- resolved production because of the strong qq component. At the LHC the 
single-resolved result is 10% of the direct while the double-resolved is reduced relative 
the single-resolved by the same factor. The LEP results suggest relatively low quark 
masses and a large resolved cross section [^. Our resolved contributions are smaller, 



possibly because the 77 luminosity at LEP drops more slowly with \J S^^ than do the 
77 luminosities for ions. 

These results are for all QQ pairs, as are the results shown for the other channels. 
However, a significant fraction of the pairs have masses below the HH threshold, espe- 
cially at RHIC. Pairs with mass M < 2mH are also produced in photo- and hadropro- 
duction. A few of these pairs will become quarkonium states ||6^. Most of them will 
hadronize into heavy-flavor hadrons, thanks to the presence of accompanying particles. 
A soft gluon can provide the energy to bring the quarks on shell. However, two-photon 
interactions occur in a vacuum with no available energy source. Pairs with M < 2mH 



may annihilate if they do not form quarkonium. Tables pCIV| and pCV| compare the total 
cross sections for all QQ pairs to those pairs with M > 2m//. In both cases, the two- 
photon cross sections are several orders of magnitude below those for hadroproduction 
and photoproduction. The cc cross sections are 0{nh) rather than mb and the bb cross 
sections are 0{ph) rather than /ib. 

The requirement M > 2mH dramatically reduces charm production. At RHIC the cc 
cross section is reduced a factor of 3 — 6 for direct 77, 4 — 16 for single resolved, and 5 — 22 
for the double-resolved. The higher LHC energy lessens the threshold effect considerably; 
the reduction is only a factor of ~ 1.9 for direct and single-resolved production and ~ 2.5 
for double-resolved production. These reductions depend strongly on the heavy quark 
mass. Our rric, 1.2 GeV, is about O.GAmD but the reduction is much smaller for larger 
charm masses. This threshold effect reduces the overall sensitivity of the calculation 
to the charm quark mass. Charm production calculations with a threshold cut match 
recent LEP data indicating that the reduced cross sections are more appropriate 



experimentally. 

The threshold reduction is smaller for bottom quarks since = 4.75 GeV ~ O.Om^. 
At RHIC, the cross section is a factor of ~ 1.5 smaller for direct photoproduction and 
1.5 — 2.5 for single- and double-resolved production. At the LHC, all the cross sections 
are reduced by (10 — 20)%. The threshold effect is more important for larger A because 
/ S^^ falls with increasing Ra- 



Figures ^ and show the ratios for QQ production with and without the threshold 
cut as a function of pt and rapidity for cc production at RHIC and bb production at 
the LHC. An invariant mass cut simply selects pairs with M > 2m//. The threshold cut 
only affects low pt quarks. The minimum pjn for a quark to pass the 2m// threshold is 



Pt > ^rn]j - m^, 1.4 GeV for m^ = 1.2 GeV and 2.3 GeV for m^ = 4.75 GeV. Naturally, 
with a larger quark mass, the minimum pt would decrease. The larger threshold effect 



for charm production appears because the peak of the px distributions in Figs. |TT] and [13 
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is below this minimum px- The average pr for bottom production is larger so that more 
of the cross section survives the threshold cut. On the other hand, the cross section is 
reduced most near yi = where low p^ dominates the rapidity distribution. At large 
rapidities, the pair QQ mass is always above threshold. The threshold has the smallest 
effect on direct production and the strongest on the double-resolved cross sections, as 
can be expected due to the decreasing effective energy available for each process. 



Energy conservation requires that a heavy quark pair produced with mass M retain 
that energy. To compensate for the 'mass deficit', M — 2m j^, the final state mesons 
must have less kinetic energy than the initial state quarks. Near threshold, the quark pt 
and rapidity distributions presented in Figs. p!T| - |T^ will differ from the final state meson 
distributions. 



We do not present any two-photon results for pA since the cross sections are very 
small and the proton substructure could play a role On the other hand, the small 
proton radius allows pA collisions to reach higher 77 energies than the corresponding AA 
collisions so that bb production would be energetically accessible in pi and pAn collisions 
at RHIC. 



The major uncertainties in the 77 calculations are the quark masses and the 77 
luminosity. In contrast to hadroproduction and photoproduction, the higher order cor- 
rections seem fairly well known. If the 77 luminosity can be sufficiently well determined, 
heavy quark production could then be used to fix the quark masses. The uncertainties 
in the 77 luminosity are comparable to those for 'yA and also depend on the minimum 
impact parameter. However, final states with known 77 couplings can be used for cal- 
ibrations. Lepton pair production covers the full range of S^j and may be accurately 
calculated using only electrodynamics. Production of well known mesons may also help 
check the luminosity. With these calibrations, (7(77 — > QQ) could clearly be measured 
to the 10% level. At that point, other theoretical uncertainties will dominate and the 
measurements can be used to determine the heavy quark masses. These masses can then 
be used in calculations of other processes. 



V. THEORETICAL COMPARISONS 



In this section, we compare and contrast some of the calculational uncertainties in 
our results. We have performed fully LO calculations, including LO parton densities and 
a one- loop evaluation of as(Q^). We first compare our LO results with NLO calculations. 
We also describe the dependence of our results on the chosen quark mass and scale. 



At any order, the partonic cross sections for all three processes studied may be 
expressed in terms of dimensionless scaling functions f-j'^^ that depend only on the 
variable rj = s/Aim? — 1 
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fc=0 i=0 




where s is the square of the partonic center of mass energy, mq is the heavy quark mass, 
and is the scale. The cross sections are expanded in powers of as(Q^) and a. The 
exponents q and p depend on the initial process: q = 1, p = 1 direct photoproduction; 
q = 0, p = 2 hadroproduction; and q = 2, p = direct two-photon production. The 
summation over k includes all corrections beyond LO which only involve powers of 
C(s{Q'^) with k = corresponding to the Born and A; = 1 to the NLO cross sections. It 
is only at NLO that the logarithms W (Q'^ / rriQ) appear. Two scales, for renormalization 
and factorization, appear in the calculation but they are generally assumed to be the 
same since this choice is used in fits of the parton densities. The total cross sections 
are obtained by convoluting the partonic cross sections with the parton distribution 
functions or photon fluxes. The scale enters the total cross section in the evaluation 
of asi^Q"^) and in the parton densities of the proton or photon (for resolved processes). 

For a fully consistent calculation of the partonic cross sections, as{Q^) should be 
evaluated to one loop when k = 0, two loops when k = 1 ■ ■ -. The strong coupling 
constant has been evaluated up to three loops, corresponding to k = 2. However, 
a consistent evaluation of the cross section, order-by-order in partonic cross section, 
parton distribution, and OsiQ"^), is usually not done. One is usually interested only in 
the effect of the next-higher-order corrections to the total cross section and it is only 
possible to measure the change by leaving other inputs, such as the parton densities, the 
same at all orders. Thus, theoretical ratios of the total NLO to LO cross sections, the 
K factors, are typically independent of the observable . 



The hadroproduction K factor is larger for 'light' heavy quarks, low pt, and low M. 
As the heavy quark mass increases, K drops from 2.5 — 3 for cc to 1.8 — 2.3 for bb in this 
energy range. For direct photoproduction, the K factors are smaller. The calculated K 
factors for direct photoproduction of bottom are 1.4 — 1.7 for yJS^p = 314 and 1265 GeV 
respectively p9[. The LO resolved photoproduction results, 0{a'j.) with p = 2, q = 
in Eq. (pi]) as at LO in hadroproduction, are typically used without NLO corrections in 
photoproduction so that the same order in asiQ"^) is used for both direct and resolved 
photoproduction. Thus the K factor would only be constant with rapidity and transverse 
momentum for direct photoproduction, not for the sum in Eq. (P). However, the resolved 
contribution is always rather small and should not significantly change the K factor. 
The NLO 77 corrections are smaller still, K = 1.33 for cc and 1.24 for bb in direct 
77, dropping to ~ 1.15 for cc and 1.21 for bb single-resolved production The 



small K factor for direct 77 should perhaps not be a surprise since, in a sense, 77 is 
the cleanest determination of the K factor because there is no parton density, as{Q^) 
or scale dependence at LO. 

This K factor, calculated with both the LO and NLO scaling functions convoluted 
with NLO parton densities and two- loop evaluations of as(Q^), does not mesh with a full 
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LO calculation using LO parton densities. The incompatibilities include the difference 
in as(Q^) evaluated at one and two loops. In the MRST LO densities, A3 = 0.204 GeV 
so that a^-^°°P = 0.287 and a^-ioops ^ q 220. The hadronic LO cross sections calculated 
with the MRST HO distributions are ctlo = 196 /ib at 200 GeV, rising to 6.03 mb at 



5.5 TeV, compared with a^LO = 382 fib at 200 GeV and 5.83 mb at 5.5 TeV. The NLO 
evaluation is two times larger at RHIC, but at the LHC, the results agree within 3%. 
The difference is mostly due to the higher one- loop as{Q'^)- Because of these variations, 
we do not apply any K factors to our LO calculations. 

Our calculations for all three processes use the same values of rriQ and Q^. The 
values are chosen to agree with hadroproduction data at fixed target energies. Pho- 
toproduction and two-photon data imply larger charm quark masses than used here. 
The typical charm mass used for those processes, 1.6 — 1.7 GeV, predict lower cross 
sections than those employing the quark masses obtained from hadroproduction. One 
can speculate as to why this is true. It may be that the incident quarks and gluons in 
hadroproduction interact with the virtual heavy quark at its current quark mass while 
the almost real photons interact with the constituent c and b quarks. The constituent 
quark mass is larger than the current quark mass and is thus more compatible with the 
photoproduction data. On the other hand, since K > 2 for hadroproduction, unincorpo- 
rated higher order corrections may explain the apparent need for different quark masses. 
Near threshold, the bb cross section has been evaluated to next-to-next-to-leading or- 
der and next-to- next-to-leading logarithm (NNLO-NNLL). Recent results from HERA-B 
69[] agree very well with the predicted 30 ± 8 ± 10 nb NNLO-NNLL cross section 



calculated with rrib = 4.75 GeV. The NLO evaluation at the same energy is a factor 
of two smaller, suggesting that NLO calculations require smaller bottom quark masses. 
However, the NNLO-NNLL expansion is invalid far away from threshold so that the 
importance of further higher order corrections is difficult to quantify. 



Figure [I^ shows the quark mass dependence of the total cross sections for all three 
processes. We plot cr(mQ)/cr(mo) where mo = 1.2 GeV for charm and 4.75 GeV for 
bottom. The scale used is oc 4m^ and ml for charm and bottom, respectively. Results 
are shown for \/S = 250 GeV Si-|-Si collisions at RHIC and 5.5 TeV Pb-|-Pb collisions 
at the LHC. The mass sensitivity is smaller at higher energies, as expected. For a 
given energy, hadroproduction is the least mass dependent. The direct photoproduction 
and two-photon production processes have nearly the same mass dependence. Resolved 
production has a stronger mass dependence, especially at RHIC. The mass dependence 
is stronger for charm than bottom, mainly because the charm mass is varied by 50%, 
from 1.2 to 1.8 GeV while the bottom mass is only varied 18%, from 4.25 to 5.00 GeV. 

The photoproduction and two-photon cross sections are more mass dependent than 
hadroproduction at the same energies. The maximum •yp collision energies, v^S'mal, 
are a factor of 4 — 6 less than v^S" in AA collisions at both colliders. The maximum 
photon energy fraction, hc/nipRA, is ~ 0.03 for gold or lead, rising to 0.1 for silicon, 
much smaller than the maximum parton energy fraction, a; = 1. This energy deficit 
is difficult to overcome and is only compensated for by the photon flux factor of 
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in charm production with the heaviest ions. The energy difference is more important 
for bottom production, especially at RHIC since y/ S^^x is close to the bb threshold. 
As Table | shows, lighter ions may be advantageous for photoproduction studies since 
the higher photon energies and accelerator luminosities can compensate for the smaller 
cross sections. The energy deficit in 77 production is even larger, 9.8 — 31.25 between 
the maximum ^JS^^ and y/S. For the heaviest ions at RHIC, the maximum yJS^^ is 

below the bb threshold. In these collisions, even charm is not far from threshold. The 
factor of Z'^ cannot compensate for such an energy deficit. The 1 mb cc cross section for 
M > 2mc at LHC is still a factor of 1000 lower than those of the other processes. Thus 
good experimental separation is essential for observing clean 77 interactions. 

The photo- and hadroproduction scale dependence is small. The cross sections only 
change a few percent between Q"^ = itiq and 4mQ because increasing the scale decreases 
OisiQ^) but increases the gluon density F^^x, Q^) and vice versa. The two effects largely 
cancel. At NLO, the scale dependence is usually larger for charm and bottom quarks 
than at LO [|r2| because ^^(Q^) multiplies the logarithm In^Q^/mQ). The scale depen- 



dence only enters through the resolved contributions in 77 production where the effect 
is a factor of 1.5 — 2 on single-resolved production and 1.05 — 1.4 on double-resolved 
production. However, the total cross section is essentially unaffected because the direct 
contribution is independent of scale. 



VI. EXPERIMENTAL SEPARATION 



To study ultra-peripheral heavy quark production, it is necessary to be able to disen- 
tangle the three channels. Photoproduction, hadroproduction, and two-photon interac- 
tions may be separated on the basis of overall event characteristics. The signatures that 
can be used to distinguish between production processes are whether there are rapidity 
gaps in the event and whether the nuclei break up. Nuclear breakup can be measured 
with downstream ZDCs. 

Other variables may be helpful in separating event classes. The event multiplicity is 
lower for photoproduction and two-photon interactions because less energy is available. 
The multiplicity also depends on the details of the interaction. For two-photon inter- 
actions, the total event px should be less than 2hc/ Ra- Unfortunately, it is necessary 
to reconstruct the entire event to measure the total pt- This is difficult in QQ events. 
Because of these difficulties, multiplicity and px will not be further considered here as a 
separation factor. 

This section will focus on isolating clean photoproduction and two-photon final states. 
The large hadroproduction cross sections are a background to these events. The cuts 
discussed here are geared toward reducing the hadroproduction background and thereby 
differentiating between the three production processes. RHIC data [0] show that almost 
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all hadronic interactions break up both nuclei. In contrast, photoproduction should 
only dissociate the target nucleus while two-photon interactions should leave both nuclei 
intact. However, the photoproduction and two-photon interactions occur at moderate 
impact parameters, less than 50 (500) fm at RHIC (LHC), where one or more additional 
photons may be exchanged. 



For heavy nuclei like lead or gold, the additional photons can excite one or both 
nuclei, leading to nuclear breakup. Except for the common impact parameter, these 
additional photons are independent of the two-photon or photonuclear events |]7D[. The 
probabilities for excitation are significant. The probability of a single given nucleus being 
excited in a collision at 6 = 2i?^ is about 30% with gold at RHIC, rising to 35% for 
lead at the LHC As 6 rises, the excitation probability drops as 1/6^. The breakup 
probability scales as A^^^^ so that for even slightly lighter nuclei like argon or silicon 
the breakup probability is considerably reduced. Since the nuclear breakup probabilities 
are independent of each other, the probability for both nuclei to dissociate is simply the 
square of the single dissociation probability. 



One way to select photoproduction events is to eliminate events where both nuclei 
break up. This should eliminate almost all of the hadropro duct ion events while retain- 
ing most of the photoproduction and two-photon interactions. Unfortunately, there are 
no calculations of the hadronic interaction cross sections without accompanying nuclear 
breakup. Indeed, such a calculation is problematic because too little is known about 
the recoil energy transfer in the nucleus. However, using a Glauber calculation, we find 
that the cross section for a single nucleon-single nucleon interaction in an Au+Au col- 
lision is about 700 mb, 10% of the total hadronic cross section. At RHIC and LHC 
energies even soft nuclear interactions involve significant energy transfer. Thus phase 
space considerations dictate that the interacting nucleons are likely to be ejected from 
the nucleus. There could then be considerable momentum transferred to the nuclear 
fragments. Here we assume that each nucleus has a 10% chance of remaining intact. 
With these assumptions, the heavy quark hadroproduction cross sections with one nu- 
cleus remaining intact are not too different from those presented in Section |T|. The 
10% chance of remaining intact may be high, but, to be conservative, we will use the 
rates in Section lll| to calculate the hadronic backgrounds to photoproduction. 



The rejection factor R is the probability of finding a rapidity gap with width y in a 
hadronic collision, R ~ exp{—ydN/dy), where dN/dy is the average multiplicity per unit 
rapidity. Here, we will only consider charged particles but with a calorimeter to detect 
neutrals, the rejection would be enhanced. The UAl collaboration parameterized the 



charged particle multiplicity as dNch/dy ~ — 0.32+ln y S/GeV Neglecting the small 
difference between pp and pp collisions, at midrapidity at RHIC, dN^h/dy ~ 2.6 



rising to dN^h/dy ~ 4.4 at the LHC. Most modern experiments use forward detectors 
like beam-beam counters to measure particle production over a wide rapidity range. 
These counters can be used to find rapidity gaps. Here we will consider two charged 
particle detectors each covering 2 units of rapidity on each side of a central detector. 
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representative of modern experiments. We scale the midrapidity multiplicities by 0.6 
because dNct/dy decreases away from midrapidity. 

For photoproduction, requiring a single rapidity gap will reject all but R = 
exp (—2 X 0.6 X 2.6) = 0.04 of the hadronic interactions at RHIC while at LHC the 
rejection for a similar detector is R = exp (—2 x 0.6 x 4.4) = 0.005. These factors, cal- 
culated for pp collisions, should be conservative for pA and AA collisions where there is 
typically more than a single nucleon-single nucleon collision. These factors would also 
apply to the rejection of photoproduction backgrounds when considering two-photon 
reactions. Since there are two rapidity gaps in two-photon interactions, these rejection 
factors are squared. 



The CDF collaboration has experimentally observed comparable rejection factors in 
a study of diffractive bottom production [|I7|. They isolated a diffractively produced bot- 
tom sample from pp collisions [|1^ with a signal to noise ratio of 3:1 despite the fact that 
diffractive production was 1/160 of the total hadronic b cross section. This corresponded 
to R = 0.002, smaller than the R values calculated above. This is probably because the 
CDF forward counters cover 2.7 units of rapidity and are supplemented with forward 
calorimeters that detect neutrals. Thus the rejection factors are not unrealistic and 
could even be improved on with more solid angle coverage. Of course, for nuclear beams 
the higher multiplicities per participant pair should increase the rejection factors, even 
for single nucleon-single nucleon interactions, due to possible secondary interactions. 



The rapidity gap requirements should reject relatively few signal events since photo- 
production always leads to a rapidity gap. The exceptions are the events that are accom- 
panied by additional electromagnetic interactions. A small fraction of these breakups 
will involve high energy photons which produce final state particles that fill in the rapid- 
ity gap, causing the signal event to be lost. Even at the LHC, the breakup probability 
due to high energy photons is quite small and events with additional particles should 
not appreciably affect the signal. If it is necessary to also require that one nucleus re- 
mains intact, then signal loss will need to be considered but such loss will not affect the 
viability of the measurement. 

Other backgrounds are neglected here. Examples include single diffractive (hadronic) 
charm production and double-Pomeron interactions. Single diffractive production will 
have one rapidity gap and could potentially mimic charm photoproduction. However, 
the diffractive final state would be produced quite near the beam rapidity, forward of 
the predominantly photoproduced charm. Double-Pomeron charm production will be in 
the central region with two rapidity gaps. Because the Pomeron has such a short range, 
both single and double diffractive interactions can only occur over a very narrow range 
of impact parameters so that their AA cross sections should be small [^. The other 
major background, beam-gas interactions, is detector and vacuum system specific and 
will not be discussed here. 



Conservatively, both the rapidity gaps and nuclear breakup criteria should reject 
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more than 99% of the hadronic events. Although these criteria are not completely 
independent, comparing the numbers in Sections || - |^ shows that application of either 
criteria should lead to a good signal to noise ratio for selecting either photoproduction 
or two-photon interactions. If 99% of the hadronic events with b > 2Ra can be rejected, 
then hadropro duct ion is only a small background to photoproduction of heavy quarks, 
one that can be controlled to the degree necessary to measure shadowing by comparing 
pA and AA. 

One could also select events with two rapidity gaps and no nuclear breakup to search 
for two-photon interactions. However, in almost all cases, the two-photon cross sections 
are a factor of at least 1,000 smaller than the photoproduction cross sections. This 
factor is smaller than the single-gap R calculated above as well as larger than rejection 
obtained by CDF and unlikely to be achieved in a real experiment. Selecting two-photon 
events may require completely reconstructing the events and using the low total pt of 
the event as a final selection criteria. This reconstruction will necessarily have a very 
low overall efficiency, thus requiring very large data sets. 

In conclusion, both rapidity gaps and the absence of nuclear breakup are effective cri- 
teria to separate photoproduction interactions. Charged-particle multiplicity and other 
event characteristics may also be useful in refining the event selection. Once a sample of 
photoproduction events have been isolated, charm and bottom production may be stud- 
ied with conventional selection techniques such as prompt leptons, separated vertices, 
and D or B meson reconstruction. 



VII. CONCLUSIONS 



We have calculated the total cross sections, quark px and y distributions, and pair 
mass spectra for hadronic, photonuclear and two-photon production of heavy quark 
pairs using a consistent set of quark masses and scales. The hadronic processes have 
the largest cross sections, followed by photoproduction and two-photon interactions. 
However, using the characteristics of rapidity gaps and nuclear breakup, photoproduction 
and two-photon processes should be cleanly separable. 

Photoproduction and two-photon measurements can be used to test the QCD calcu- 
lations. Shadowing has a 10% effect on cc production with heavy nuclei at RHIC, rising 
to 20% at the LHC By comparing the production cross sections in pA and AA collisions, 
most theoretical uncertainties cancel so that shadowing can be accurately measured if the 
photon flux is well known. By using appropriate calibration signals, it appears that the 
photon flux uncertainties can be understood and useful gluon shadowing measurements 
made. 
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TABLES 







A A 

AA 


pA 




A 


VSnn (GeV) 


Caa (cm ^s 


VSnn (GeV) 


LpA (cm ^s 






RHIC 









250 


9.8 X 10^^ 


250 


1.2 X 10^° 


Si 


250 


4.4 X 10^® 


250 


8 X 10^^ 


I 


208 


2.7 X lO^^^ 


208 


2 X 1029 


Au 


200 


2 X 10^6 


200 


6 X 10^8 






LHC 









7000 


1.6 X 10^9 


9900 


1.0 X 10^^ 


Ar 


6300 


4.3 X 10^^ 


9390 


5.8 X 10^° 


Pb 


5500 


4.2 X 10^6 


8800 


7.4 X 1029 


TABLE L Luminosities and beam energies for AA and pA collisions at RHIC and LHC. 






AA 






A 


^beam (GeV) 


IL femax (GeV) 


^max (TeV) 


V-^max (GeV) 






RHIC 









125 


133 12.7 


2.31 


66 


Si 


125 


133 8.5 


1.92 


60 


I 


104 


111 3.9 


0.81 


39 


Au 


100 


106 3.2 


0.66 


35 






LHC 









3500 


3730 357 


1820 


1850 


Ar 


3150 


3360 185 


1080 


1430 


Pb 


2750 


2930 87.0 


480 


950 






pA LHC 









4950 


5270 435 


3630 


2610 


Ar 


4700 


5000 276 


2410 


2130 


Pb 


4400 


4690 139 


1220 


1500 


TABLE II. Beam energies, -Ebeamj Lorentz factors, 7^, photon cutoff energ 


;y in the center of 


mass 


frame, fcmax, and in 


the nuclear rest frame, £^max, 


and equivalent nucleon-nucleon center 


of mass energy, V'S'max, for AA collisions at RHIC and the LHC. Since the ion and proton beam 


energ 


;ies are expected to be the same in pA and AA collisions at RHIC, we 


only distinguish 



the pA energies at LHC. 
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AA 


(j° (d* = 1) (mbj 


cr'^''^(EKS98) (mb) 


a^^yb^ = 1) (mbj 


cr ''^(EKS98j (mbj 






RHIC 






0+0 


0.067 


0.068 


0.0019 


0.0019 


Si+Si 


0.30 


0.31 


0.0080 


0.0083 


I+I 


8.96 


9.74 


0.199 


0.206 


Au+Au 


24.8 


27.4 


0.530 


0.550 






LHC 






0+0 


2.35 


2.13 


0.351 


0.346 


Ar+Ar 


23.3 


20.4 


3.00 


2.93 


Pb+Pb 


1790 


1500 


190.0 


186.7 


TABLE III. Direct and resolved cc photoproduction cross sections integrated over b > 2Ra 


in peripheral AA collisions at RHIC and LHC. 






AA 




a'i''^(EKS98) (//b) 




CT-^^^ (EKS98) {nh) 






RHIC 






0+0 


0.047 


0.049 


0.0031 


0.0031 


Si+Si 


0.178 


0.188 


0.0116 


0.0115 


I+I 


2.33 


2.46 


0.154 


0.148 


Au+Au 


4.94 


5.22 


0.332 


0.315 






LHC 






0+0 


11.9 


11.4 


2.93 


2.93 


Ar+Ar 


107 


102 


22.2 


22.6 


Pb+Pb 


718 


686 


121 


126 


TABLE IV. Direct and resolved bb photoproduction cross sections integrated over b > 2Ra 


in peripheral AA collisions at RHIC and LHC. 








cc 




bb 




pA 


^dir^^i = 1) (^b) 






^res(^i = 1) (nb) 






RHIC 






pO 


2.68 


0.081 


2.34 


0.154 


pSi 


7.29 


0.213 


5.79 


0.378 


pi 


54.2 


1.32 


23.5 


1.52 


pAu 


100 


2.34 


35.5 


2.31 






LHC 






pO 


110 


20.7 


630 


202 


pAr 


485 


83.9 


2670 


774 


pPh 


7940 


1190 


40100 


9910 



TABLE V. Direct and resolved cc and bb photoproduction cross sections integrated over 
b> rp + Ra in pA collisions at RHIC and LHC. 
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cc 



bb 



AA 


N{S' = 1) 


iV(EKS98) 


N{S' = 1) 


iV(EKS98) 






RHIC 






0+0 


6.75 X 10^ 


6.94 X 10*5 


4.88 X 10^ 


5.10 X 10^ 


Si+Si 


1.36 X lO'^ 


1.41 X 10"^ 


8.35 X 10^ 


8.73 X 10^ 


I+I 


2.47 X lO'^ 


2.69 X lO'^ 


6.70 X 10^ 


7.06 X 10^ 


Au+Au 


5.07 X 10^ 


5.60 X 10^ 


1.06 X 10^ 


1.10 X 10^ 






LHC 






0+0 


4.15 X 10*^ 


3.80 X IQS 


2.29 X 10^ 


2.20 X 10^ 


Ar+Ar 


1.13 X 10^ 


9.98 X 10^ 


5.58 X 10^ 


5.37 X 10^ 


Pb+Pb 


8.29 X 10^ 


7.05 X 10^ 


3.58 X 10^ 


3.46 X 10^ 



TABLE VI. Total cc and bb photoproduction rates in peripheral AA collisions over a 10 s 



run at RHIC and LHC. The rates are based on Tables III and IV 





cc 


bb 


pA 


N{S' = 1) 


N{S' = 1) 




RHIC 




pO 


3.30 X 10^ 


2.99 X 10^ 


pSi 


6.00 X 10^ 


4.93 X 10^ 


pi 


1.11 X 10^ 


5.00 X 10^ 


pAu 


6.08 X 10^ 


2.27 X 10^ 




LHC 




pO 


1.32 X 10^ 


8.32 X 10^ 


pAi 


3.11 X 10^ 


1.88 X lO'^ 


pPh 


6.70 X 10^ 


3.69 X lO'^ 



TABLE VII. Total cc and bb photoproduction rates in pA collisions over a 10^ s run at 
RHIC and LHC. The rates are based on Table 0. 



AA a{S' = 1) (mb) cj(EKS98) (mb) (t(EKS986) (mb) 

RHIC 

0+0 4.04 3.93 4.00 

Si+Si 8.54 8.22 8.42 

I+I 22.6 21.6 22.3 

Au+Au 33.1 31.6 32.6 

LHC 

0+0 113 93.2 104 

Ar+Ar 426 323 379 

Pb+Pb 1090 714 948 

TABLE VIII. Total cc hadroproduction cross sections integrated over b > 2Ra in periph- 
eral A A collisions at RHIC and LHC. 
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AA 



N{S' = 1) 



iV(EKS98) 



Af(EKS986) 



0+0 
Si+Si 
I+I 
Au+Au 

0+0 
Ar+Ar 
Pb+Pb 



3.96 X 10^ 
3.76 X 10^ 
6.10 X 10^ 
6.62 X 10^ 

1.74 X IQio 
1.83 X 10^° 
4.57 X 10^ 



RHIC 



LHC 



3.85 X 10^ 
3.62 X 10^ 
5.84 X 10^ 
6.33 X 10^ 

1.44 X 10^° 
1.39 X 10^° 
3.00 X 10^ 



3.92 X 10^ 
3.70 X 10^ 
6.02 X lO'^ 
6.52 X 10^ 

1.61 X 10^° 
1.63 X 10^° 
3.98 X 10^ 



TABLE IX. Total cc hadroproduction rates in peripheral AA collisions at RHIC and LHC 



with a 10^ s run. The rates are based on Table IVIHI. 



AA 



aiS^ = 1) (^b) 



(t(EKS98) (^b) 



fT(EKS98fe) (^b) 



0+0 
Si+Si 
I+I 
Au+Au 

0+0 
Ar+Ar 
Pb+Pb 



22.7 
47.9 
111 
154 

2490 
9110 
21700 



RHIC 



LHC 



24.1 
51.7 
127 
183 

2260 
7930 
17500 



23.3 
49.5 
117 
167 

2390 
8600 
20200 



TABLE X. Total bb hadroproduction cross sections integrated over b > 2Ra in peripheral 
AA collisions. 



AA iV(y = 1) iV(EKS98) Af(EKS986) 

RHIC 

0+0 2.22 X 10^ 2.36 x 10^ 2.27 x 10^ 

Si+Si 2.11 X 10^ 2.27 x 10^ 2.18 x 10^ 

I+I 3.01 X 10^ 3.43 X 10^ 3.16 x 10^ 

Au+Au 3.17 X 10^ 3.67 x 10^ 3.34 x 10^ 

LHC 

0+0 3.84 X 10^ 3.48 x 10^ 3.68 x 10*^ 

Ar+Ar 3.92 x 10^ 3.41 x 10^ 3.70 x 10*^ 

Pb+Pb 9.11 X 10*^ 7.35 x 10^ 8.48 x 10^ 

TABLE XI. Total bb hadroproduction rates in peripheral AA collisions with a 10^ s run at 
RHIC and LHC. The rates are based on Table 
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A 

pA 


a[b = 1) (mb) 


cc 

(7(EKS98) (mb) 


00 

aiS' = 1) (^b) 


(7(EKS98) (/ib) 






RHIC 






pO 


4.19 


4.14 


23.7 


24.4 


pSi 


7.33 


7.21 


41.5 


43.1 


pi 


26.2 


25.8 


131 


139 


pAu 


38.6 


38.1 


187 


201 






LHC 






pO 


153 


138 


3740 


3540 


pAi 


368 


318 


8850 


8180 


pPh 


1820 


1460 


43000 


38200 



TABLE XII. Total cc and bb hadroproduction cross sections in minimum bias (all b) pA 



collisions at RHIC and LHC. 



cc 



bb 



pA 


N{S' = 1) 


Af(EKS98) 


N{S' = 1) 


Af(EKS98) 






RHIC 






pO 


5.02 X 10^ 


4.96 X 10^ 


2.85 X 10"^ 


2.92 X lO'^ 


pSi 


5.86 X 10^ 


5.77 X 10^ 


3.32 X lO'^ 


3.45 X lO'^ 


pi 


5.25 X 10^ 


5.17 X 10^ 


2.62 X 10"^ 


2.79 X lO'^ 


pAu 


2.32 X 10^ 


2.28 X 10^ 


1.12 X lO'^ 


1.20 X lO'^ 






LHC 






pO 


1.53 X 10^2 


1.39 X 10^2 


3.75 X 10^° 


3.57 X 1010 


pAi 


2.01 X 10^2 


1.75 X 10^2 


4.87 X 10^° 


4.50 X 10^° 


pPh 


1.35 X 10^2 


1.07 X 10^2 


3.18 X 10^° 


2.82 X 10^° 



TABLE XIII. Total cc and 66 hadroproduction rates in minimum bias pA collisions over a 



10 s run at RHIC and LHC. The rates are based on Table XH 
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AA 


a (nb) 


n — 

ah cc 
a (nb) 


a (nb) 


0" (nb) 


M > 2m£) 
a (nb) 


a (nb) 








RHIC 








0+0 


4.64 


0.08 


0.0020 


1.65 


0.022 


0.00039 


Si+Si 


32.0 


0.49 


0.013 


10.8 


0.125 


0.0023 


I+I 


1320 


10.7 


0.345 


288 


1.18 


0.027 


Au+Au 


3650 


22.2 


0.786 


601 


1.37 


0.035 








LHC 








0+0 


236 


11.7 


0.24 


128 


6.01 


0.10 


Ar+Ar 


4530 


210.0 


4.36 


2410 


105 


1.76 


Pb+Pb 


1110000 


45000 


951 


565000 


21400 


352 


TABLE XIV. Two photon cc cross sections in peripheral AA cohisions at RHIC and LHC, 



integrated over b > 2Ra- 



all bb M > 2mB 

AA a'^" (pb) a^-''''' (pb) a^-'""^ (pb) a'^" (pb) fj^"''^" (pb) a'^'''''' (pb) 

RHIC 

0+0 0.268 0.0018 0.00038 0.194 0.0010 0.00029 

Si+Si 0.923 0.0031 0.00083 0.582 0.0013 0.00046 

LHC 

0+0 285 31.7 3.08 262.6 28.9 2.62 

Ar+Ar 4890 491.0 49.3 4480 444 41.7 

Pb+Pb 943000 75000 8260 855000 66800 6820 

TABLE XV. Two photon bb cross sections in peripheral A A collisions at RHIC and LHC, 
integrated over b > 2Ra- 



all QQ 



M > 2mH 



AA 


N{cc) 


N{bb) 


N{cc) 


N{bb) 






RHIC 






0+0 


4.62 X 10^ 


2.65 X 10^2 


1.64 X 10^ 


1.91 X 10^2 


Si+Si 


1.43 X 10^ 


4.08 X 10^2 


4.81 X 10^ 


2.57 X 10^2 


I+I 


3.61 X 10^ 




7.81 X 10^ 




Au+Au 


7.36 X 102 


LHC 


1.21 X 10^ 




0+0 


3.97 X 10^ 


5.11 X 10^ 


2.15 X 10^ 


4.71 X 10^ 


Ar+Ar 


2.04 X 10^ 


2.33 X 10^ 


1.08 X 10^ 


2.13 X 10^ 


Pb+Pb 


4.84 X 10^ 


4.41 X 10^ 


2.47 X 10^ 


3.90 X 10^ 



TABLE XVI. Total cc and bb two-photon rates in peripheral AA collisions over a 10 s run 



at RHIC and LHC. The rates are based on Tables XIV and XV 
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FIG. 2. Charm photoproduction in peripheral Au+Au colUsions at RHIC for b > 2Ra. 
The single c quark pT (a) and rapidity (b) distributions are shown along with the cc pair 
invariant mass (c) . The direct (dashed) , resolved (dot-dashed) , and the sum of the two (solid) 
are shown. The direct contribution is divided by two to distinguish it from the total while 
the resolved contribution is multiplied by ten. There are two curves for each contribution: 
S"* = 1 and EKS98. At this energy, the curves arc almost indistinguishable but the curves with 
shadowing are somewhat higher, especially at negative rapidities. In the rapidity distributions, 
the photon is coming from the left. 
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FIG. 3. Bottom photoproduction in peripheral Au+Au collisions at RHIC for b > 2Ra- 
The single b quark pT (a) and rapidity (b) distributions are shown along with the bb pair 
invariant mass (c). The direct (dashed), resolved (dot-dashed), and the sum of the two (solid) 
are shown. The direct contribution is divided by two to distinguish it from the total. There 
are two curves for each contribution: 5' = 1 and EKS98. At this energy, the curves are almost 
indistinguishable. In the rapidity distributions, the photon is coming from the left. 
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FIG. 4. Charm photoproduction in peripheral Pb+Pb collisions at LHC for b > 2i?^. The 
single c quark pT (a) and rapidity (b) distributions are shown along with the cc pair invariant 
mass (c). The direct (dashed), resolved (dot-dashed), and the sum of the two (solid) are 
shown. The direct contribution is divided by two to distinguish it from the total. There are 
two curves for each contribution: 5* = 1 and EKS98. The unshadowed curves are higher 
than the shadowed, particularly at large rapidities. In the rapidity distributions, the photon 
is coming from the left. 
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FIG. 5. Bottom photoproduction in peripheral Pb+Pb collisions at LHC for b > 2Ra- 
The single b quark pT (a) and rapidity (b) distributions are shown along with the bb pair 
invariant mass (c). The direct (dashed), resolved (dot-dashed), and the sum of the two (solid) 
are shown. The direct contribution is divided by two to distinguish it from the total. There 
are two curves for each contribution: = 1 and EKS98. The unshadowed curves are higher 
than the shadowed. In the rapidity distributions, the photon is coming from the left. 
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FIG. 6. Charm hadroproduction in peripheral Au+Au coUisions at RHIC for b > 2Ra- The 
single c quark pT (a) and rapidity (b) distributions are shown along with the cc pair invariant 
mass (c). The curves are 5* = 1 (solid), EKS98 (dashed), and EKS986 (dot-dashed). At this 
energy, the three curves are almost indistinguishable. 
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FIG. 7. Bottom hadroproduction in peripheral Au+Au collisions at RHIC for b > 2i?yi. 
The single b quark (a) and rapidity (b) distributions are shown along with the bb pair 
invariant mass (c). The curves are 5* = 1 (solid), EKS98 (dashed), and EKS986 (dot-dashed). 
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FIG. 8. Charm hadroproduction in peripheral Pb+Pb collisions at LHC for h > 2Ra- The 
single c quark px (a) and rapidity (b) distributions are shown along with the cc pair invariant 
mass (c). The curves are 5* = 1 (solid), EKS98 (dashed), and EKS986 (dot-dashed). 
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FIG. 9. Bottom hadroproduction in peripheral Pb+Pb collisions at LHC for b > 2i?^. The 
single b quark pT (a) and rapidity (b) distributions are shown along with the bb pair invariant 
mass (c). The curves are 5* = 1 (solid), EKS98 (dashed), and EKS986 (dot-dashed). 
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FIG. 11. Charm production by two-photon processes in peripheral Si+Si colhsions at 
RHIC. The results are shown for all pairs with no mass cut. The single c quark pT (a) 
and rapidity (b) distributions are shown along with the cc pair invariant mass (c). The solid 
curve is the sum of all contributions: a'^^'^ (dashed), a^~^^^ (dot-dashed), and a^"^^^ (dotted). 
The direct contribution is divided by two to facilitate comparison. Since either photon can 
be resolved, the single-resolved rapidity distribution is reflected around yi = to account for 
both sources. 
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FIG. 12. Bottom production by two-photon processes in peripheral Si+Si collisions at 
RHIC. The results are shown for all pairs with no mass cut. The single b quark pT (a) and 
rapidity (b) distributions are shown along with the bb pair invariant mass (c) . The solid curve 
is the sum of all contributions: (j^'"^ (dashed), a^~'^^^ (dot-dashed), and a'^~'^^^ (dotted). The 
direct contribution is divided by two to facilitate comparison. Since either photon can be 
resolved, the single-resolved rapidity distribution is reflected around yi = to account for 
both sources. 
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FIG. 13. Charm production by two-photon processes in peripheral Pb+Pb collisions at 
LHC. The results are shown for all pairs with no mass cut. The single c quark pT (a) and 
rapidity (b) distributions are shown along with the cc pair invariant mass (c) . The solid curve 
is the sum of all contributions: a"^^"^ (dashed), a^~'^^ (dot-dashed), and a^""^^ (dotted). The 
direct contribution is divided by two to facilitate comparison. Since either photon can be 
resolved, the single-resolved rapidity distribution is reflected around yi = to account for 
both sources. 
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FIG. 14. Bottom production by two-photon processes in peripheral Pb+Pb coUisions at 
LHC. The results are shown for all pairs with no mass cut. The single b quark pT (a) and 
rapidity (b) distributions are shown along with the bb pair invariant mass (c) . The solid curve 
is the sum of all contributions: (j^'"^ (dashed), a^~'^^^ (dot-dashed), and a'^~'^^^ (dotted). The 
direct contribution is divided by two to facilitate comparison. Since either photon can be 
resolved, the single-resolved rapidity distribution is reflected around yi = to account for 
both sources. 
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FIG. 15. Reduction in charm production due to the requirement that M > 2mD in 
two-photon production in peripheral Si+Si coUisions at RHIC. The ratio of the cross section 
above threshold relative to the total cross section is shown as a function of pT (a) and yi (b) 
for cr^'"' (sohd), a^''^^^ (dashed) and a^~^^^ (dot-dashed). Since either photon can be resolved, 
the single-resolved rapidity ratio is reflected around yi = to account for both sources. 



50 




1.0 — 




— 1.0 



— 0.8 




(b) 



— 0.6 



0.0 







1 2 



3 



-5.0 



-2.5 



0.0 



2.5 



5.0 



Pt (GeV) 



yi 



FIG. 16. Reduction in bottom production due to the requirement that M > 2mB in 
two-photon production in peripheral Pb+Pb colhsions at LHC. The ratio of the cross section 
above threshold relative to the total cross section is shown as a function of px (a) and yi (b) 
for a^^^ (sohd), a^~'^'^^ (dashed) and cr^"'^^® (dot-dashed). Since either photon can be resolved, 
the single-resolved rapidity ratio is reflected around yi = to account for both sources. 
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FIG. 17. The quark-mass dependence of our calculated cross sections. The left side is for 
charm, normalized to the cross sections with rric = 1.2 GeV while the right side is for bottom, 
normalized to cross sections with = 4.75 GeV. The RHIC results are for Si+Si interactions 
while the LHC results are for Pb+Pb interactions. The photoproduction ratios in (a) and 
(b) are for direct (solid-RHIC; dot-dashed-LHC) and resolved (dashed-RHIC; dotted-LHC) 
production. The hadroproduction results at RHIC and LHC are given by the solid and dashed 
curves respectively in (c) and (d). The two-photon ratios in (e) and (f) are for a'^^^ (solid-RHIC; 
dotted-LHC), a^''^ (dashed-RHIC; dot-dot-dot-dashed-LHC) and a^-'^ (dot-dashed-RHIC; 
dot-dash-dash-dashed-LHC) . 
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